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Abstract. The Faddeev-Yakubovsky equations in configuration space have 
been solved for the four nucleon system with special interest in their scat- 
tering states. We present results concerning the structure of the first 4 He, scat- 
tering states in different (T,S) channels including the first inelatic threshold 
N+3N^NN+NN and the n-t cross section using realistic potentials. 

1 Introduction 

The theoretical description of the A=4 continuum constitutes, in addition to 
its technical difficulty, a serious challenge for the NN interaction models. The 
remarkable success encountered in the description of the 3N system and in the 
4N bound state Q is largely based on the three nucleon interaction (TNI) . The 
two body dynamics alone looses « 10% in the 3N binding energy and » 15% 
in the 4N. TNI compensate this defect but, a part from an energy scaling, are 
found to have small influence in the low energy scattering states Contrary 
to the A=3 case, these states present several nearthreshold structures, which 
can be hardly sensible to the TNI and should be reproduced with the two-body 
forces only. A first example is provided by the n+t cross section which manifests 
a much more vivid behaviour than the flat n-d one, with a resonant structure 
at Tiab w 3 MeV (figure |l|). A second interesting accident is the existence of 
quasidegenerate n- 3 He and p- 3 H thresholds with, in their middle, the first 4 He 
excitation and all that in an energy gap of 0.76 MeV. The position of such a 
state, which dominates the inter-threshold dynamics, has to be ensured with 
an accuracy better than 100 KeV in order to reproduce the experimental cross 
sections, what is a strong requirement for an excited state. 

2 Formalism 



In the case of four identical particles interacting via a two-body potential V, the 
Faddeev-Yakubovsky (FY) equations result into two integrodifferential equa- 
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Figure 1. Comparison between the nd and nt elastic cross sections (in barns) 



tions, coupling two amplitudes denoted K and H: 

(E-H -V)K = V [(P 23 + Pis) (e + P34) K + e(P 23 + Pis) H] 
(E-H -V)H = V [(P 13 P 2i + P14P23) K + P13P24 H] (1) 

in which Py are the permutation operators and e = ±1 depending on we deal 
with bosons or fermions. The total wavefunction is reconstructed according to: 

SP = $1+3 + ^2+2 

f 1+ 3 - [l + e(Pi3 + P 23 )] [1 + e(Pi4 + P 24 + P34)] K (2) 
^2+2 = [1 + e(Pi 3 + P23 + Pi4 + P 24 ) + P3P24] H (3) 

Each amplitude F=K,H considered as a function of its own set of Jacobi coor- 
dinates, is expanded in the corresponding angular variables according to 

(xyz\F)=Yl ^^Y a {x,y,z). (4) 

F a are the unknowns and Y a are tripolar harmonics containing spin, isospin and 
angular momentum variables. Label a holds for the set of quantum numbers 
defined in a given coupling scheme and includes the type of amplitude K or H. 

Boundary conditions for scattering states are of Dirichlet-type. In the 1+3 
elastic case e.g. they are implemented by imposing at large enough value of z 

K(x,y,z) = t(x,y) 
H(x,y,z) = 

t(x,y) being the triton Faddeev amplitudes previously determined. FY equa- 
tions ensure a solution which, for n+t S-waves, behaves asymptotically like 

K(x,y,z) - t(x, y) sin (qz + S) 

where S is the phaseshift and q is related to the center of mass kinetic energy 
T cm and to the physical momentum k by T cm = ^q 2 = f-^-fc 2 . A more detailed 
explanation of the formalism and numerical methods used can be found in |j. 
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3 Results 

We have been interested in the first + excitation of 4 He, experimentally man- 
ifested as a r—0.5 MeV resonant state, 0.4 MeV above the pt threshold Q]. 
Using MT I-III potential || we got j|] the binding energies (B) and r.m.s. 
radius (R) displayed in table [j], what means a bound state with B<i ffe , =0.26 



Table 1. Binding energy (MeV) 
and r.m.s. radius (fm) 

4 He 4 He* 3 H 



B 
R 



30.30 8.79 
1.44 4.95 



8.53 
1.72 



Table 2. Low energy N+3N parameters (fm) 
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MeV below the N+NNN threshold. Using AV14 || and Nim II potentials 
we found close results and similar conclusions were also reached using sepa- 
rable interactions To investigate the structure of this state, we looked at 
the two-body correlation functions. The results (figure ||) show a superposition 
of two structures with different length scales, the short distance part being 
similar to the triton one. This suggests - rather than a breathing mode as it 
is widely considered in the literature - a 3+1 structure formed by an almost 
unperturbed triton plus a fourth nucleon orbiting around. This state, with a 
r.m.s. R « 5 fm, constitutes actually the first " halo- nucleus" . It seems difficult 
for a strong interaction to generate a first excitation in the continuum and the 
open question is whether or not Coulomb forces, when added to the existing 
models, would be able to put it in the right place with the required accuracy. 

Within the same model we have calculated (table |2|) the parameters of the 
effective range expansion 



;(fc) = k cot{5) 
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- + \r k 2 + v Q k 4 + o(fc 6 ) 
a 2 



(5) 



for N+NNN S-states. Expansion (||) provides an accurate representation of 
the below breakup scattering although the existence of the 4 He* pole makes 
necessary the explicit inclusion of a singular term for the T=S=0 channel. 
It is worth noticing the coherence between the He* binding energy and the 
scattering parameters from table [|. If expansion (||) is inserted in the S-wave 
scattering amplitude /o _1 (^) = g(k) — ik, 2 imaginary poles follow and the 
nearest to threshold k- = 0.095 fm" 1 corresponds to E- = —0.25 MeV, in 
agreement with the direct calculations of the 4 He* binding. 

The 4N scattering states with two open channels N+NNN^NN+NN has 
been calculated for the S=T=0 final state ||. Of particular interest is the 
extraction of the imaginary part of the strong NN+NN scattering length which 
controls the fusion rate in the process d+d — > n+ 3 He. We found qr = +4.91 ± 
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Figure 2. Two-body correlations for (a) triton, (b) 4 He ground and (c) first excited 
states 4 He*. Solid (dashed) line denotes the triplet (singlet) contributions. In (d), the 
results of 4 He* are compared to the triton correlation function suitably scaled. 



0.02 fm and a very small value of aj = — 0.0115±0.0001 fm which should be only 
slightly modified once the Coulomb interaction is switched on. This value is due 
to the small overlapping between the K and H configurations which respectively 
govern the N+NNN and NN+NN asymptotic states. Other calculations of the 
dd cross section at very low energy exist || but no values of a/ were given. 

Special attention was paid to the n+t scattering. We found |3| that the 
description provided by the simple MT I-III model fits remarkably well the 
data, as it can be seen in figure || where the differential cross sections at different 
energies are displayed. Using realistic potentials, we have performed fic|l the 
same calculations for the J v = + , 1 + , 0~, 1~, 2~ states with 1 S'o, 3 Si, 3 D\ 
waves in Vnn and expansion (^) limited to l x , l y , l z = 0, 1, 2. The singlet (ao), 
triplet (ai), coherent (a c ) scattering length and zero energy cross section (ao) 
are in table [| compared with experimental values taken from jnj. They are 
in very good agreement with those obtained by Pisa group fl2]| ao = 4.32 and 
ai = 3.80 but the infered a c and cr(0) values are well far from experiment. 



Table 3. n+t scattering length with realistic interactions 

ao ai a c = jap + |ai c(0) = 7r(ap + Sa\) 

AV14 4~31~ 3J9 3J2 lfM 

Nijmll 4.31 3.76 3.90 192 

MT I-III 4.10 3.63 3.75 177 

exp 3.59±0.02 170± 3 
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Calculations have been pursued beyond zero energy. One can see from figure 
^ (dashed curve), that two-body realistic potentials alone fail in describing the 
low energy cross section as they failed in giving the binding energies B3 and B4. 
We would like to emphasize that this failure is in fact already present in the n+d 
case, although somehow hidden because the sizeable disagreement concerns 
only the doublet scattering length which turns to be very small compared to 
the quartet one and it is furthermore suppressed by statistical factors. 
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Figure 3. n-t differential cross section with Figure 4. n-t total cross section with 
MT I-III at different n laboratory energies AV14 compared to experiment 



To improve this situation we introduced the phenomenological TNI: 

w( P ) = w r w a — p = V% 2 + y 2 

p p 

Using the parameter set W r = 500 W a — 1 74 ji = 2.0, we got an overall low 
energy agreement with the values #3=8.48 MeV, #4=29.0 MeV, ao=4.0 fm, 
ai = 3.53 fm, a c — 3.65 fm and cr(0)=168 fm 2 and a satisfactory n+t S-wave 
cross section (solid line). We point out that some dispersion exists in the exper- 
imental results which makes difficult the comparison with theoretical values. If 
one assumes the two last measurements [[ll] of a c to be valid, our calculations 
are not far from the proposed value a c = 3.60. In this case, however, our scat- 
tering lengths are not compatible with those infered from experiment, which 
are in their turn incompatible with each other. A more precise measurement of 
er(0) would be helpful. For that purpose the possibilities offered by the neutron 
beam facility, recently approved at CERN jl3]], could deserve some attention 
from the nuclear Few-Body community. 

If the inclusion of TNI is determinant in the n+t S-wave region, it doesn't 
improve the resonance peak, on the contrary. In a first tentative to explain this 
discrepancy we have included P- waves in the NN interaction pO). We found 
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that although separated P-waves had significant influence in the cross section, 
their global effect was very small and they could not explain the disagreement. 
A recent work done by A. Fonseca H] using AGS equations and 1-rank approx- 
imation, was able to push expansion (||) further than we could do and found a 
satisfactory description. It would be interesting to have independent confirma- 
tion of this result and go beyond the first obstacle of the A=4 continuum. 

Acknowledgement. The numerical calculations have been performed with the 
Cray-T3E of the CGCV (CEA) and IDPJS (CNRS). 
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